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Section  I 


INTRODUCTION 


This  report  presents  the  results  of  a program  undertaken  with  the 
Naval  Air  Development  Center  (NADC)  to  analyze  flight  test  data  obtained 
from  a Navy  YT-2B  aircraft,  instrumented  and  operated  by  NADC.  The  objectives 
of  this  program  are  to  estimate  the  longitudinal  stability  and  control  param- 
eters of  the  YT-2B,  from  digitally  recorded  flight  data,  collected  and  supplied 
by  NADC,  using  an  iterated  Kalman  filter/ fixed-point  smoother  parameter  identi- 
fication algorithm.  A linear,  constant  coefficient,  small  perturbation  equations- 
of-motion  model  is  used  to  represent  the  characteristics  of  the  YT-2B. 

The  data  consists  of  a total  of  seventeen  flight  records:  five 

records  at  a low  speed  landing  approach  configuration  and  twelve  at  a higher 
speed  configuration  with  Mach  number  around  .6.  From  this  group,  three  low 
speed  and  five  high  speed  records  were  initially  selected  as  the  mos'*  promising 
for  further  data  analysis.  The  results  reported  herein  are  primarily  those 
obtained  using  the  low  speed  flight  records.  These  results  correlated  favorably 
with  each  other  and  with  the  initial  parameter  estimates  supplied  by  N'ADC. 
Inaccuracies  in  the  estimated  parameters  and  differences  between  parameters 
extracted  from  different  flight  records  are  discussed  and  straightforward 
solutions  posed. 


The  organization  of  this  report  is  as  follows:  A brief  description 
of  the  identification  algorithm  and  the  mathematical  model  of  the  YT-2B  is 
first  presented,  followed  by  some  general  comments  related  to  the  flight  data. 
The  identification  results  and  conclusions  and  recommendations  follow.  Tables 
and  figures  are  presented  at  the  end  of  the  report  for  the  reader's  convenience. 
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Section  II 


IDENTIFICATION  ALGORITHM 


The  general  problem  of  estimating  stability  and  control  parameters 
from  flight  test  records  can  be  viewed  as  one  of  fixed-point  smoothing,  that 
is,  estimating  both  the  unknown  parameters  and  the  initial  condition  of  the 
aircraft  from  the  data.  With  this  as  a starting  point,  the  estimation 
philosophy  may  then  be  divided  into  two  major  groups:  Bayesian  and  non-Bayesian. 
The  difference  between  these  groups  is  that  Bayesian  philosophy  can  take  into 
account  a priori  information  about  the  data  of  unknowns,  whereas  non-Bayesian 
does  not.  The  Bayesian  philosophy  assumes  that  the  entire  information  available 
to  an  estimator  is  contained  in  the  a posteriori  density  function;  that  is, 
the  density  function  of  the  unknowns  given  the  data. 


side red: 


With  the  Bayesian  philosophy,  two  estimation  criteria  can  be  con- 

(1)  The  minimum  variance  or  minimum  mean  squared  error  criterion; 
the  resulting  estimator  is  the  conditional  mean,  that  is,  the 
expected  value  of  the  a posteriori  density  function,  or 

(2)  Maximum  likelihood  (Bayesian)  or  most  probable  estimate,  which 
is  the  mode  of  the  a posteriori  density  function. 


The  identification  algorithm  used  here  was  developed  in  Reference  1 
using  the  mean  squared  error  criterion  with  appropriate  Gaussian  assumptions. 

By  annexing  a constant  parameter  vector  model  of  unknowns  to  the  aircraft  state 
equations,  the  recursive  algorithm  which  results  is  a nonlinear  iterated 
Kalman  filter/fixed-point  smoother.  The  nonlinear  filter  is  a form  of  an  ex- 
tended Kalman  filter  utilizing  a "local  iteration"  scheme  between  successive 
measurements  to  reduce  errors  in  linearizing  the  equations  of  motion  and  mea- 
surement system  models  by  improving  the  reference  trajectory.  This  improve- 
ment is  accomplished  by  smoothing  each  measurement  data  point  backwards  in 
time  one  point  and  rc- 1 inearizing.  In  addition,  the  outputs  of  the  filter  at 


each  data  point  can  he  used  in  a fixed-point  smoothing  algorithm  to  produce  a 
better  smoothed  estimate  of  the  unknown  parameters  and  initial  condition  of  the 
aircraft.  The  algorithm  is  summarized  in  Appendix  A. 

This  technique  is  both  a response  error  and  an  cqua* ion  error  mini- 
mization technique  under  the  usual  Gaussian  assumptions  in  the  error  sources. 
That  is,  the  technique  adjusts  the  parameters  of  the  model  so  as  to  minimize 
the  "weighted"  errors  between  both  the  measured  responses  (accelerations 
included)  and  the  error  in  formulating  the  equations  of  motion  of  the  airplane. 
The  weighting  in  the  measurements  is  selected  to  be  compatible  with  the 
accuracy  of  the  instrumentation  and  recording  system,  whereas  the  weighting 
of  the  equation  error  is  compatible  with  the  error  in  the  mathematical  model 
of  the  airplane.  These  error  sources  are  commonly  referred  to  as  the  meas- 
urement and  process  noise,  respectively.  Although  these  sources  of  error 
are  modeled  as  white,  Gaussian  and  independent  in  the  identification 
algorithm  with  covariances  R and  Q , icrpcctivcly , they  can  also  be  con- 
sidered as  artificial  approximations  to  non-random  type  errors  such  as  non- 
linearities  in  the  model  or  lags  in  the  instrumentation  system,  for  example, 
and  as  such  serve  as  effective  weighting  matrices. 

The  selection  of  these  noise  statistics  in  particular  is  an  iterative 
process  which  is  carried  out  in  practice  by  observing  the  residual  sequences 
(measurement  data  minus  predicted  measurements)  of  the  filtering  operation  and 
adjusting  the  statistics  as  required  to  obtain  consistency  between  the  pre- 
dicted and  actual  dispersion  of  the  residuals.  If  the  model  is  correct 
(implies  both  form  and  statistics  used),  the  actual  dispersion  of  the  residuals 
should  be  zero  mean,  white  and  consistent  with  their  theoretically  calculated 
statistics.  Two  indicators  arc  available: 

(1)  Visual  inspection  of  the  residuals  in  the  form  of  time 
history  plots;  and 

(2)  Estimation  of  f?  and  Q from  the  residuals  using 

A A 

covariance  matching  techniques  (£  and  Q). 
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t?  and  Q are  estimated  by  the  covariance  matching  technique  by  approximating 
the  actual  residual  covariance  by  its  "sample  covariance"  and  by  using  the 
relationships  between  this  residual  covariance  and  the  ft  and  Q matrices 
in  the  identification  algorithm.  If  the  model  is  correct  and  the  records 
sufficiently  long,  these  estimates  serve  as  a form  of  hypothesis  test  to 
determine  if  the  initial  R and  Q were  selected  properly  or  if  not,  a guide 
to  further  "tuning." 

Two  other  pieces  of  information  are  required  to  perform  the  identi- 
fication. These  are: 

• Initial  estimates  of  parameters  and  states 

• Variance  of  the  initial  estimate  ( P0  ). 

The  selection  of  these  items  are  discussed  in  the  results  section. 
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Section  III 


MATHEMATICAL  MODEL  AND  FLIGHT  RECORDS  ANALYZED 


Model  of  YT-2B 

The  rigid  body,  three-degree-of- freedom  linearized  longitudinal 
equations  of  Motion  model  used  to  describe  the  recorded  small  perturbation 
■otions  of  the  YT-2B  for  identification  purposes  is  given  in  Table  1.  Also 
included  are  the  Measurement  system  model  and  appropriate  definitions.  This 
model  is  essentially  equivalent  to  the  one  supplied  by  NADC  (which  is  the 
standard  "BuAer  type"  stability  axes  systems  model)  except  that  the  units  are 
different  and  the  unknown  stability  and  control  parameters  to  be  identified 
are  completely  aerodynamic  in  origin.  In  particular,  the  units  for  angular 
measurements  are  degrees,  linear  velocities  and  linear  accelerations  have  the 
units  of  feet  and  the  and  v stability  parameters  do  not  include 

kinematic  terms.  The  moment  parameters  are  primed,  for  example  , indi- 
cating that  the  effects  of  Cm.  have  been  lumped  into  these  parameters. 

, 

V , M 0 and  Za  are  used  to  compensate  for  out  of  trim  or  improperly  chosen 
reference  conditions. 

As  shown  in  Table  1,  eight  measurements  are  modeled.  These  consist 
of  perturbation  variables  in  pitch  attitude  ( true  airspeed  ( A\/m  ) , 

pitch  rate  ( ),  angle  of  attack  ( ),  ^ -stability  axis  linear  ac- 
celeration ( ),  two  angular  acceleration  measurements  ( q i , ) 

7 m ' m ' * m 

and  v -stability  axis  linear  acceleration  Bias  parameters, 

which  are  subscripted  with  a small  b,  are  also  included  in  each  measurement. 


The  perturbation  measurements  supplied  by  NADC  were  formed  by  sub- 
tracting an  appropriate  reference  condition  from  each  total  flight  measurement 
and  using  the  following  transformation  to  obtain  the  perturbation  stability  axis 
linear  acceleration  measurements: 


c.osa^ 


- sm  o' 


cosec. 


a~  (t») 


(i0) 
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where  nv  , arc  the  body  reference  axes  accelerometer  measurements,  is 
trim  angle  of  attack  and  ),  ( t0  ) are  used  to  develop  perturbation 

variables. 

Flight  Records  Analyzed 


Seventeen  flight  records  were  supplied  by  NADC.  For  reference  pur- 
poses, these  consisted  of  five  low  speed  records  in  the  landing  approach  con- 
figuration, numbered  2 through  6,  and  twelve  at  a higher  speed  configuration, 
labeled  7 through  19  (except  15).  Of  these,  three  low  speed  records  (2,  5, 
and  6)  and  five  high  speed  records  (10,  12,  13,  18,  and  19)  were  initially 
selected  as  the  most  promising  for  identification  purposes.  The  flight  con- 
ditions for  each  record  are  summarized  in  Table  2.  All  records  were  at  a 
sample  rate  of  20  samples  per  second. 

The  initial  values  for  the  coefficients  of  the  equations  of  motion 
were  supplied  at  Vt  « 236  fps  and  Vt  - 679  fps.  These  initial  values  were 
converted  to  the  parameters  of  the  model  given  in  Table  1 and  also  updated 
to  be  closer  to  the  actual  flight  condition  for  comparison  purposes.  These 
values,  which  served  as  initial  estimates,  are  given  in  Table  3.  In  all  cases, 
'Z'f  was  assumed  equal  to  zero. 
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Section  IV 


GENERAL  COMMENTS 

Prior  to  discussing  the  results  in  the  next  section,  a few  general 
comments  about  the  flight  data  and  a discussion  of  problem  areas  which  were 
uncovered  in  analyzing  the  low  speed  records  are  in  order.  'Ibis  will  provide 
a proper  perspective  for  interpreting  the  identification  results.  Specifically, 
four  potential  problem  areas  are  worthy  of  discussion  since  these  are  con- 
sidered partly  responsible  for  inaccuracies  in  the  estimated  parameters  and 
the  differences  between  parameters  exti acted  from  different  flight  records. 

The  first  is  the  uncertainty  of  the  actual  elevator  control  input 
applied  to  the  airplane.  It  is  well  known  that  the  control  input  time  history 
recorded  for  these  flight  records  is  extremely  noisy  and  is  not  the  control 
input  applied.  This  is  readily  apparent  when  observing  that  the  flight 
variables,  such  as  q , a and  ^ , for  example,  in  most  of  the  records  have  for 
all  practical  purposes  reached  a steady  state  but  the  measured  control  input 
has  not.  The  difficulty  is  not  so  much  the  high  frequency  content  of  the  con- 
trol input  noise  but  low  frequency  components  which  are  within  the  bandwidth 
of  the  airplane  and  have  a signal  strength  which  is  sufficiently  large  compared 
to  the  actual  control  input.  These  low  frequency  components  can  be  observed 
by  comparing  the  flight  responses  generated  using  the  identified  parameters 
in  the  model  with  the  flight  data,  by  filtering  the  control  input  to  reduce 
the  high  frequency  noise  components  or  by  observing  the  residuals  in  the  iden- 
tification technique.  If  is  not  filtered,  the  high  frequency  components 
do  of  course  affect  the  measurement  noise  levels  which  should  be  used  on  the 
acceleration  measurements,  particularly  £ ^ and  . This  error,  of 
course,  is  a modeling  error  which  should  properly  be  construed  as  a source  of 
process  noise. 


For  example,  the  measured  control  input,  » can  be  modeled  as 

the  true  control,  be  , plus  additional  noise,  ur  , as  shown  in  Equation  (2) 
below.  Solving  Equation  (2)  for  Se 


1 


se  - or 


(2) 


* 


'm 


and  substituting  the  results  in  the  pitching  moment  equation  yields 


9 = * MSe  Sem  + MSe  uT  (3) 

process  noise 

where  the  M<,£uf  term  represents  the  equation  error  or  process  noise  due  to 
the  error  in  measuring  the  control  input.  The  next  problem  is  to  determine 
if  the  noise  term  is  white  or  colored  and  its  variance.  If  it  is  colored,  then 
a shaping  filter  driven  by  white  noise  can  be  used  to  model  the  error  term 
where  the  parameters  of  the  filters  are  themselves  unknown  constants  to  be 
identified  or  else  they  could  be  modeled  as  random.  Assuming  that  the  noise 
is  stationary  and  its  characteristics  do  not  change  with  the  application  of 
a control  input,  the  characteristics  of  the  noise  term  can  be  obtained  from 
the  sanplc  power  spectral  density  (PSD)  of  the  5e/y>with  and  without  the  control 
input  applied  if  the  record  is  sufficiently  long  to  provide  meaningful  results. 
For  example,  this  was  done  for  records  6 and  13,  but  the  complete  time  history 
was  used  for  records  2,  5 and  18.  The  results  are  given  in  Figures  la  through 
lg.  These  PSD  plots  serve  as  a guide  in  the  selection  of  the  process  noise 
level  to  account  for  the  noise  in  the  control  input  measurement. 


It  should  be  realized  that  or  is  acting  like  an  unknown  forcing 
input  to  the  airplane  model,  similar  to  an  unknown  gust  input.  Consequently, 
in  generating  time  history  responses  to  compare  with  the  flight  data,  an  esti- 

A 

mate  of  the  unknown  forcing  function  time  history,  ur  ( t ) , should  be  included 
in  the  model,  especially  if  uf  ( t ) has  very  low  frequency  components.  An 
estimate  of  ur{t)  is  directly  available  in  the  identification  technique  if 

A 

a shaping  filter  is  used  in  the  model,  as  described  above.  If  not,  ur  (t) 
can  be  generated  using  the  algorithm  developed  in  Reference  1. 

Of  course,  the  objectives  here  are  to  identify  the  stability  and  con 
trol  parameters  of  the  YT-211,  not  the  characteristics  of  the  noise  corrupting 
the  control  input  time  history  unless  it  becomes  necessary  to  do  so.  Less 
sophisticated  approaches  include  the  following: 
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(1) 


-r  J— i-. Sr  -tt-s  - 


Do  nothing.  Use  the  Se  as  measured  and  adjust  the  measure- 
m'nt  noise  statistics  accordingly. 


(2)  Similar  to  (1),  except  in  addition  use  a constant  level  of 
white  process  noise  as  indicated  by  the  PSD  plots  of  the 
control  input. 

(3)  Similar  to  (1),  except  first  filter  Se  to  eliminate  the 
high  frequency  components  of  noise. 


(4)  Similar  to  (3)  above,  but  in  addition  use  a constant  level  of 
white  process  noise  to  condensate  for  the  low  frequency  content 
of  the  noise.  The  level  is  initially  established  from  the 
PSD  plot  of  the  control  input  measurements  as  described  above. 


(5)  Conjecture  about  what  the  actual  control  input  was  and  modify 
accordingly. 


Because  of  the  limited  scope  of  this  investigation  and  other  possible 
instrumentation  difficulties  discussed  below,  approaches  3 and  5 were  taken 
initially,  followed  by  approach  2.  In  particular,  the  control  input  on  the  low 
speed  records  was  filtered  with  a "nonrealizable"  digital  filter  which  has  the 
frequency  characteristics  of  two  cascade  fourth-order  Butterworth  filters,  but 
with  no  phase  shift,  with  a cutoff  frequency  of  3.5  Hz.  For  record  6,  it  was 
also  conjectured  that  the  control  was  actually  equal  to  a constant,  before  and 
after  the  application  of  the  pulse  type  input,  and  was  therefore  set  equal  to 
the  sample  mean  of  the  input.  The  actual  and  filtered  inputs  are  shown  in  Fig- 
ure 2 for  comparison.  Also  shown  are  the  inputs  filtered  with  a cutoff  fre- 
quency of  2 Hz.  A common  cutoff  frequency  of  3.5  Hz  was  selected  by  observing 

the  PSD  plots  of  the  control  inputs  so  as  not  to  destroy  the  signal  content  and 

also  to  preserve  the  sharp  edges  of  the  pulse  type-  input  of  record  6.  Clearly, 

a lower  frequency  could  be  used  for  records  2 and  5.  In  actuality,  little 

differences  were  noted  in  the  "major"  parameters  identified  for  the  low  speed 
records  with  the  use  of  the  filtered  or  unfiltered  control  input,  if  the 
measurement  noise  statistics  were  adjusted  to  compensate  (improperly)  for  this 
added  error. 
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For  the  high  speed  records,  this  problem  is  more  serious  because  of 
the  smaller  control  input  amplitudes  with  the  same  level  of  noise  and  the  in- 
creased bandwidth  of  the  airplane. 

The  second  problem  is  concerned  with  an  inconsistency  between  the 
normal  stability  axis  linear  acceleration  measurement  ( Aa^.)  and  the  angle 
of  attack  vane  measurement  ( a-v  ).  This  inconsistency  is  very  apparent  when 
examining  the  residuals  in  the  identification  technique  or  by  comparing 
generated  responses  with  the  flight  data.  It  has  the  features  of  a phase  shift 
between  the  two  measurements  or  a possible  phase  lag  and  scale  factor  error 
in  the  vane  measurement.  The  error  is  definitely  nonrandom  and  the  re- 
sulting parameter  estimates,  specificaliy  , Mg g and  XA  , 

are  dependent  upon  the  noise  levels  assigned  to  these  measurements.  This 
presents  difficulties  in  establishing  the  proper  weighting  for  these  measure- 
ments (measurement  noise)  because  they  must  now  compensate  for  the  non-random 
nature  of  thase  errors.  However,  one  would  suspect  the  error  to  be  in  the  a.y 
measurement  and  indications  are  given  below  to  support  this  suspicion, 
although  conclusive  evidence  is  lacking. 


A third  problem  is  associated  with  the  use  of  stability  axes  pertur- 
bation variables,  especially  the  fla*  stability  * -axis  linear  acceleration 
measurement.  It  turns  out  that  the  derivation  of  an  accurate  measure- 

ment from  the  Y)^  and  Yt^  body  reference  axes  linear  accelerometers  requires 
an  accurate  knowledge  of  trim  angle  of  attack  ( cct  ).  The  transformation 
is  given  by  Equation  (1).  Usually,  during  the  short  period  portion  of  the 
flight  record,  / A j >>  / j , and  from  Equation  (1)  it  is  seen  that  an 
error  in  <ti  can  cause  large  errors  in  Aa^  during  this  portion  of  the  record. 

This  problem  was  detected  in  an  attempt  to  match  the  Aa.^  measure- 
ment through  the  short  period  portion  of  the  records.  Fair  matches  were 
obtained  with  the  inclusion  of  an  Xg  parameter,  but  the  value  extracted 
was  unrealistically  large.  For  example,  the  results  given  in  the  following 
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table  were  obtained.  The  value  of  Xy  identified  was  around  the  same  order 


Record  2 Records  5 and  6 


Para. 

Initial 

Estimated 

Initial 

5 

Estimated 

6 

Estimated 

X*  -.205 

-.386 

-.172 

-.25 

-.27 

*$  0 

.325 

0 

.14 

.55 

of  magnitude  as  X# 

, which  ic  absurd! 

It  is  now 

clear  that 

the  Se  time  his 

tory  looked  sufficiently  close  to  the  error  introduced  by  the  transformation 
(that  is,  the  s in  eet  Y>*  term)  to  partially  condensate  for  its  effect. 
Excluding  the  Xg£  parameter  and  using  the  measurement,  could,  of 

course,  introduce  errors  in  the  X- force  parameters  extracted,  depending  upon 
the  level  of  measurement  noise  statistics.  Adjusting  the  measurement  noise 
assigned  to  A a*.  to  condensate  for  this  error  can  be  done,  but  if  a large 
value  is  required  with  respect  to  the  other  measurement  noise  statistics, 
not  using  the  measurement  at  all  provides  the  same  effect.  Therefore,  the 
latter  has  been  chosen.  Unfortunately,  this  reduces  the  ability  to  extract  the 
Xy  and  stability  derivatives  accurately. 

A possible  error  in  is  also  supported  by  two  additional  pieces 

of  information.  At  the  beginning  of  each  flight  record,  where  the  airplane  was 
known  to  be  near  trim  and  level  flight  with  7"  ■ 0,  the  pitch  attitude  and  abso 
lute  angle  of  attack  were  differing  by  as  much  as  3.5  degrees.  Since 

T - B - a.  -for  <f>  - 0, 

an  offset  in  0 , cc  or  both,  is  indicated.  Additionally,  from  information 
supplied  by  NADC  (Reference  2)  concerning  the  instrumentation  calibration  it 
was  noted  that  the  inclination  of  the  YT-2B  nose  boom  was  calibrated  with 
respect  to  u "mean  wing  chord  line"  and  not  the  fuselage  reference  line. 
Although  the  angular  difference  is  small,  this  could  be  a source  of  error  if 
the  linear  accelerometer  package  is  aligned  to  the  fuselage  reference  line. 
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The  difficulties  in  resolving  the  aforementioned  problem  are  for  the 
*0!  t part  due  to  the  use  of  "perturbation"  flight  variables,  which  were  supplied 
and  hence  used  in  the  equations  of  motion  model  and  measurement  system  model. 

From  Table  1,  it  is  seen  that  in  actuality  there  is  no  mechanism  to  identify 
a true  offset  in  the  ocy  measurement  so  as  to  determine  oct  more  accurately. 

The  use  of  a different  model  is  therefore  recommended  for  future  flight  test 
data  analysis  at  these  flight  conditions.  One  such  model  which  would  suffice 
employs  total  flight  variables,  has  linear  dimensional  aerodynamics  (including 
thrust)  but  nonlinear  kinematics  and  incorporates  a nonlinear  time  varying 
stability  to  body  axis  transformation  in  the  measurement  system  to  model  the  body 
reference  axis  accelerometer  measurements  ( r>x  and  ) directly.  With  ad- 
ditional biases  modeled  in  the  measurement  system  to  be  identified  as  unknowns, 
this  model  will  provide  a mechanism  by  which  to  identify  an  offset  in  CcY  so  as 
to  rectify  the  above  problem.  For  example,  this  model  is  briefly  outlined  in 
Equations  (4)  and  (5)  below. 


Equations  of  Motion: 

V =■  Xa  - g sin  2T 
oi>  = Za  + ^ cos  7f  + (f. 


Acceleration  Measurements: 


10  x ~ X, a.  Co5  * ~ Vs  in  %l  « Xa-  Vac 


= X.  sinoc  *■  7.  V 


» X oi  + 2 V 

A A 


Xa  - X0  ♦ Xv  AV  * X#  Aoc  * *tc. 


= ZY  AY  + * etc. 


Ma  - + M^AV  * etc. 


AV  - V - V0  , etc. 


7 = G-oc 
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The  last  area  is  an  observed  phase  lag  in  the  angular  accelerometer 
with  the  application  of  sharp  control  inputs  like  those  in  record  6.  This  is 
in  direct  agreement  with  the  large  phase  shift  observed  and  reported  in  Refer- 
ence 2 when  calibrating  this  instrument  at  NADC  on  the  oscillating  rate  table. 
The  specification  for  this  instrument  indicates  a natural  frequency  in  the 
order  of  10  Hz,  which  is  low  for  an  acceleration  sensor,  especially  when 
used  to  record  the  motions  due  to  sharp  control  inputs. 

In  general,  all  four  problems  are  more  serious  for  the  high  speed 
flight  records  because  of  the  increased  bandwidth  of  the  airplane  at  this  flight 
condition.  The  results  are  discussed  next. 


Section  V 


t 


4 


RESULTS 

This  section  presents  the  identification  results  obtained  using 
records  2,  5,  and  6.  Additionally,  the  results  from  n rather  limited  analysis 
of  records  13,  18  and  19  are  given.  Also  cited  are  the  problems  discussed  in 
the  preceding  section.  In  all  cases,  unless  otherwise  indicated,  the  initial 
parameter  estimates  are  those  given  in  Table  3.  Parameters  which  have  very 
small  effects  on  the  flight  responses,  and  therefore  cannot  be  identified  accu- 
rately, were  held  fixed  during  the  identification. 

For  example,  and  were  set  equal  to  zero  - was 

only  identified  during  assessment  of  the  problems  associated  with  the  Ad# 
measurement,  as  discussed  previously.  Also,  it  was  assumed  that  the  airplane 
was  trimmed  at  7"t  * 0,  so  that  and  were  known  exactly.  Even  if 

4 0,  the  latter  two  parameters  are  extremely  small  and  have  no  apparent 
effect  on  the  flight  responses,  especially  for  the  short  period  type  records. 
Similarly,  for  a rigid  airplane,  dCLl 9V  = 0,  so  that  x?y  depends  only  upon 
the  knowledge  of  Vf  . Since  its  value  is  somewhat  insensitive  to  small 
variations  in  Vt  , it  also  was  held  fixed  at  the  particular  estimated  trim 
condition. 


For  all  cases,  biases  ir  pitch  rate  ( ) and  all  the  acceleration 

measurements  were  identified  along  with  the  out-of-trim  parameters,  Y-0  , Mg 
and  Z0  . Since  the  measurements  used  are  perturbation  variables  for  a selected 
reference,  these  biases  were  extracted  differently  for  each  record  (as  antici- 
pated), and  therefore  are  not  presented.  The  first  data  points  in  each  record 
were  used  as  the  initial  state  estimates  (tempered  with  engineering  judgment 
because  of  the  high  level  of  noise)  and  the  variances  of  these  estimates  were 
set  equal  to  the  measurement  noise  variance  of  the  respective  state  measure- 
ments (biases  included).  The  variances  of  the  initial  parameter  estimator 
( "PQ  ),  which  represent  the  accuracy  of  these  estimates  in  the  form  of  "a 
priori"  information,  were  set  sufficiently  large  to  insure  that  the  initial 
estimate  had  little  effect  on  the  final  extracted  values. 
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The  measurement  noise  statistics  were  initially  established  by  means 
of  a "sample  variance"  calculation  on  each  measurement  for  all  records.  The 
mean  value  of  each  measurement,  which  is  required  for  such  a calculation,  was 
"approximated"  by  using  the  same  digital  filter  as  was  applied  to  the  control 
input.  The  results,  averaged  overall  records,  with  anomalies  eliminated,  are 
given  in  Table  4.  These  values  provide  an  "approximate"  estimate  to  the  random 
conponent  of  measurement  error  and  therefore  serve  as  a good  initial  estimate, 
more  like  a lower  bound,  since  they  do  not  account  for  such  things  as  instrument 
non linearities  or  dynamics. 

It  was  noted  that  the  noise  on  the  angle-of-attack  vane  for  record 
2 was  higher  than  records  5 and  6,  which  is  consistent  with  the  pilot  comments 
indicating  the  presence  of  moderate  turbulence  for  this  record.  Assuming  the 
absence  of  turbulence  in  records  5 and  6,  the  standard  deviation  of  the  verti- 
cal component  of  the  gust  velocity,  <7^  , in  the  body  axes  coordinate  system, 

was  estimated  to  be  approximately  equal  to  .5  feet  per  second.  This  value  could 
serve  as  a good  initial  estimate  in  the  selection  of  the  process  noise  level 
in  using  a turbulence  model,  such  as  a simplified  Dryden  model,  for  example, 
to  account  for  the  effects  of  the  turbulence  input  during  the  identification. 
However,  the  effects  of  turbulence  at  this  level  on  the  identification  results, 
if  they  are  neglected,  are  probably  small  compared  to  the  other  problems  cited 
above  and,  therefore,  a turbulence  model  was  not  employed. 

Although  unknown  forcing  inputs  to  the  airplane  or  model  of  the  air- 
plane are  known  to  exist  because  of  the  presence  of  turbulence  and  the  uncer- 
tainty of  the  control  input  applied,  many  of  the  results  reported  herein  assume 
no  process  noise  ( Q = 0).  Additionally,  the  angular  acceleration  measurement 
( A^)  which  was  derived  from  the  tail  and  nose  linear  accelerometers,  was  not 
used  in  the  identification,  but  the  time  history  matches  to  this  measurement  are 
presented. 
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Record  Number  2 (Q~0) : 


* 


The  identification  results  for  this  record  are  given  in  Table  5 and 
typical  time  history  response  comparisons  with  the  flight  data,  using  the  esti- 
mated parameters  for  the  particular  cases  indicated,  arc  given  in  Figures  3 
and  4.  In  all  cases,  the  control  input  was  filtered  using  a cutoff  frequency 
of  3.S  Hz  as  explained  previously.  This  record  has  little  phugoid  excitation 
(small  velocity  change)  so  that  parameters  such  as  XY  and  Af'  would  be  hard  to 
identify. 

Table  5 presents  the  identification  results  of  five  cases,  labeled 
2-1  through  2-5,  and  the  measurement  noise  statistics  used  to  perform  the  iden- 
tification. These  five  results  serve  to  highlight  some  of  the  problems  dis- 
cussed in  the  preceding  section.  The  first  case  (2-1)  uses  the  AaK  measure- 
ments with  identified,  whereas  the  other  four  cases  do  not,  to  show  the 

i improve  me  nt  gained  in  matching  the  response  in  the  former  case.  Unfor- 

tunately, is  difficult  to  identify  accurately  without  the  use  of  this  mea- 
surement . 


For  the  second  case  (2-2),  the  <S  -values  obtained  from  the  "sample 
variance"  calculation  (Table  4) , were  used  for  the  measurement  noise  covariance 
matrix  ( /?  ) , with  X g held  fixed  at  zero.  The  response  comparisons  with 
the  flight  data  and  the  residuals  (measurement  data  minus  the  predicted  esti- 
mates) during  the  identification  are  shown  in  Figures  3a  and  3b  respectively 
for  this  case.  The  unknown  forcing  input  to  the  model  (caused  mainly  by  the 
uncertainty  of  the  control  input),  the  inconsistency  of  the  Aoc  and  Aa^  mea- 
surements and  the  problems  associated  with  A<z^  are  readily  apparent.  The 
residuals  show  a definite  breakup  or  deterministic  type  dispersion,  a clear 
indication  of  a modeling  problem.  Figure  4a  gives  the  response  comparison  with 
for  case  2-1  where  is  identified.  Comparing  this  with  Figure  3a 

shows  an  indication  of  the  problem  with  A<XX  . Although  not  shown,  the  use 
of  the  Ad*  measurement  without  identifying  Xg ■ provides  a similar  response 

match  to  Aay  as  shown  in  Figure  3a. 


In 


Cases  2-3  and  2-4  show  the  effects  of  adjusting  the  values  of  the 
measurement  noise  statistics,  particularly  2-4,  where  the  accuracy  of  Act  is 
decreased  and  that  of  AO-.^  increased.  Changes  in  some  of  the  major  parameters 
arc  evident.  Figure  4b  gives  the  response  comparisons  with  Ait  and  A(l^  for 
case  2-4;  the  other  responses  were  similar  to  case  2-2.  Comparison  of  Figure 
4b  with  the  Act  and  Aa ^ response  matches  in  Figure  3a  provides  another  indi- 
cation of  the  Aec  and  Aa^  inconsistency.  Pickup  of  large  wind  components 
by  the  angle-of-attack  vane  was  not  considered  in  this  comparison  because  the 
results  were  similar  for  record  5,  which  was  obtained  in  only  slight  turbulence 
on  a different  day. 

For  the  last  case  (2-S),  the  measurement  noise  statistics  were 
selected  to  be  equal  to  the  diagonal  elements  from  the  estimator  for  the  mea- 
surement noise  statistics  ( f?  ) to  make  the  two  consistent.  This  is  not 
exactly  proper,  since  R reflects  errors  introduced  by  the  unknown  input,  which 
is  a process  noise  error,  and  other  non-randem  type  instrument  errors  and  there- 
fore, the  estimated  values  are  expected  to  be  different  for  each  record.  How- 
ever, it  provides  a convenient  guide  for  the  selection  of  R . 

The  accuracy  of  the  parameter  estimates,  as  given  by  twice  the  square 
root  of  the  diagonal  elements  of  the  final  covariance  ( 2 crf  ) , are  presented 
in  the  last  column  of  Table  5.  Case  2-5  results  were  used,  because  the  value 
of  the  measurement  noise  statistics  is  more  representative  of  the  errors  pre- 
sent in  the  flight  record.  Although  these  values  are  only  approximations  be- 
cause of  errors  which  were  modeled  inproperly,  they  appear  to  be  very  reasonable 
indications  of  the  estimation  accuracy,  which  is  construed  as  ±Z(X^. 

Record  Number  5 iO  - O) : 

Similar  to  the  identification  results  presented  for  record  2,  a 
spectrum  of  five  identification  cases  is  presented  for  record  5.  The  only  dif- 
ference is  that  Xv  was  held  fixed  at  the  initial  estimate  because  of  the  small 
velocity  change  in  th  s record.  The  results  arc  given  in  Table  6,  for  the  five 
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cases  labeled  5-1  through  5-5,  which  arc  in  direct  correspondence  to  eases  2-1 
through  2-5. 

The  response  comparisons  with  the  flight  data  and  the  residuals 

A 

are  given  for  cases  5-2  ( k!  -sample  variance)  and  5-5  ( ki  ) in  figures  5 and  7, 
respectively.  For  comparison  witli  these  time  history  matches.  Figure  6a  gives 
the  response  conparison  wi.n  Aa.y  for  case  5-1,  and  Figures  6a  and  6b  gives  the 
response  comparisons  with  Aoc  and  Aa^  for  cases  5-3  and  5-4,  respectively. 
Conclusions  drawn  from  these  results  are  similar  to  those  briefly  discussed 
for  record  2 above. 

Record  Number  6 (Q  - 0)  • 

This  record,  unlike  the  former  two,  has  a classical  pulse  type 
input  for  excitation  of  the  phugoid.  Again  five  identification  cases  are  pre- 
sented in  Table  7,  labeled  6-1  through  6-5  and  various  response  comparisons  to 
the  flight  data  are  given  in  Figures  8 through  11.  In  all  cases,  Acl%  was  not 
used. 

The  first  two  cases  (6-1  and  6-2)  used  only  the  first  five  second 
short  period  portion  of  the  flight  record  and  therefore,  Xy  was  held  fixed. 
Also,  the  control  input  was  filtered  with  a cutoff  frequency  of  3.5  Hz.  Case 
6-1  employed  the  sample  variance  computation  for  R whereas  R was  adjusted 
upwards  for  6-2.  The  response  comparisons  with  the  flight  data  are  shown  in 
Figure  8 for  case  6-2  and  indicate  the  phase  lag  in  discussed  above,  which 
is  the  main  reason  for  presenting  these  cases.  Also,  the  error  in  the  control 
input  is  clearly  indicated  in  the  q and  £ comparisons. 

Cases  6-3  and  6-4  use  a ’nodified  control  input  and  the  first  14 
seconds  of  the  flight  record.  It  was  conjectured  that  the  elevator  surface 
was  actually  a constant  before  and  after  the  application  of  the  pulse,  so  it  was 
set  equal  to  its  average  value  at  these  times.  There  is  a danger  in  doing  this 
because  of  the  presence  of  an  actual  input,  so  that  only  the  first  14  seconds  of 
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this  record  was  used.  It  appeared  as  if  an  input  were  applied  near  this  time. 
Case  6-3  employs  the  sample  variance  calculation  for  & , with  the  exception 

9 

that  04  was  adjusted  upward  to  account  for  the  phase  lag  in  . Case  6-4  uses 
the  diagonal  element  from  the  R (measurement  noise  covariance  estimator). 

There  is  no  guarantee,  however,  that  this  was  the  control  input  applied.  The 
response  comparison  with  the  flight  data,  which  is  similar  for  both  cases,  and 
the  residuals,  are  given  in  Figure  9. 


The  flight  responses  for  case  6-3  were  then  compared  for  the  first 
28  seconds  using  a step  size  {At  ) of  0.1  seconds.  Drift  was  noted  in  A&  , 
An.  and  AV  . Consequently,  Vy,  and  My  were  identified  using  the  full  28 
second  portion  of  the  record  with  the  other  parameters  held  fixed  at  the  case 
6-3  estimator  with  At  » .1.  The  results  are  presented  as  case  6-5  in  Table  7 
and  the  response  comparisons  with  the  flight  data  are  given  in  Figure  10.  A 
drift  is  still  noted  in  the  AcL  response  which  could  be  the  result  of  an  ap- 
plied input  or  a minor  parameter  which  was  not  estimated,  such  as  or  , 
for  example.  Figure  11  gives  the  response  conparison  with  the  flight  data 
using  the  "actual"  control  input  for  case  6-5  but  at  a At  = .05  seconds  for 
comparison  purposes. 


Record  Number  5(Q/0) 


Record  5 was  reanalyzed  using  approach  number  2 described  previously. 
That  is,  process  noise,  w,  was  used  to  model  the  uncertainty  of  the  control 
input  as  described  in  Equation  (3).  The  process  noise  was  assumed  zero  mean, 
white  and  Gaussian  with  a variance  level  selected  from  the  PSD  plot  of  the  con- 
trol input,  Figure  lb.  It  was  assumed  that  the  high  frequency  portion  of  the 
PSD  represented  noise  only,  the  characteristics  of  which  were  similar  at  the 
low  frequency  end.  Th>'  dashed  line  in  Figure  lb,  which  is  representative  of 
the  level  selected,  gives  a process  noise  standard  derivation  of 


. deg 


where  the  factor  of  2 is  required  because  of  the  definition  used  between  the 
autocorrelation  function  and  the  power  spectral  density  in  the  PSD  computer 
program 
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The  identification  results  employing  this  level  of  process  noise, 
an  unfiltered  control  input  and  a level  of  measurement  noise  consistent  with 
the  R estimates  arc  given  in  Table  8 as  case  5-6.  Only  the  short  period  and 
control  effectiveness  parameters  were  identified;  the  others  were  held  fixed 
at  the  initial  estimates.  The  results  compare  favorably  with  those  for  Q = 0 
and  6e  filtered  (Table  6)  and  the  ± 2 cr^  values  appear  to  be  very  reasonable 
indications  of  the  estimation  accuracy. 


For  Case  5-6,  the  time  histories  of  the  state  estimates  and  the 
residuals  are  given  in  Figure  12.  Since  an  estimate  of  the  unknown  forcing 
function,  &■(*),  was  not  available  for  incorporation  into  the  model  for  genera- 
ting time  history  responses  to  compare  with  the  flight  data,  none  is  given. 
However,  the  state  estimates  from  the  iterated  Kalman  filter  and  the  residuals 
suffice  because  they  include  the  effects  of  process  noise.  Similar  to  the 
results  given  above,  the  Zla  and  Aa residuals  still  show  an  inconsistency. 
Additionally,  the  white  noise  approximation  for  may  be  questionable,  although 
more  analysis  is  required  to  determine  this,  partly  because  of  the  cause  and 
effect  relationships  due  to  the  Aac  and  Aau„  uncertainty.  If  a shaping  filter 
is  used  to  color  the  sample  rate  of  20  per  second  may  be  low. 


Record  Number  5 {Aac.  and 


inconsistency) : 


As  cited  above  and  described  previously,  an  inconsistency  appeared 
to  exist  between  the  Aac  and  the  measurement  and  this  error  made  it  diffi- 

cult to  establish  a consistent  set  of  measurement  noise  statistics.  From  analysis 
on  record  5,  it  appeared  like  a scale  factor  error  in  the  vane  measurement  could 
be  responsible.  Usually,  such  an  error  is  the  result  of  wing/fuselage  (possibly 
landing  gear)  air  flow  effects  and  is  largest  for  short  booms  and  low  Mach 
numbers.  Unfortunately,  the  measurement  system  model  in  Tabic  1 did  not  include 
a scale  factor  parameter  which  could  be  treated  as  an  unknown  to  be  identified. 


\ 
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However,  a model  was  available  from  previous  identificati  >n  efforts 
at  Calspan  (Reference  3)  which  uses  the  kinematic  equations  of  the  aircraft  and 
a typical  measurement  system  model  in  the  manner  of  a strapped-down  inertial 
measuring  unit.  This  model,  when  used  in  the  identification  algorithm,  allowed 
for  the  estimation  of  such  parameters  as  biases  and  scale  factors  in  the  instru- 
mentation system,  "independent"  of  the  aerodynamic  model  of  the  airplane.  This 
is  desirable  here  because  of  the  modeling  error  caused  by  the  noisy  control  input. 
The  model  is  given  in  Equations  (6)  and  (7)  below. 

Dynamic  Model : 


9 

s 

9m  - n * n 

0 

m 

u 

s 

- gsin  0 b An*  — tig 

■ 

au.  a cos  0 + An,  - n. 
r 7 )b 

t *3 

Measurement  System: 

9m 

K 

9 + +V4 

*em 

S 

0 + Ofi  + 

£ 

U * 

ab  + ^ 

S 

(oc2-f  Mrz)t/2  -/■  vb  b ir7 

These  are  the  three -degree-of- freedom  kinematic  equations  written  in  the  air- 
craft body  axis  coordinate  system,  where  the  subscripts  m , b , and  s imply 
measurements,  biases,  and  scale  factors,  respectively.  Noise  terms  are  repre- 
sented by  ir,  through  v'f  . 

For  application  to  record  5,  biases  were  used  to  relate  the  total 

flight  variables  to  the  perturbation  measurements.  Additionally,  An,  and 

An,  were  approximated  by  Aa,  and  Aa~  , respectively,  and  the  noise 
iw  n?  i m 

statistics  used  are  those  given  in  Table  4,  with  the  exception  of  . OZ 

was  increased  to  account  for  the  above  approximation.  Five  oiases,  ^ , 

^ib  ’ 9 b an<*  °^b  am*  sca*e  ^actor»  Ks  * were  treated  as  unknowns.  The 
result  was  a Ks  estimated  value  of  .275--surprisingly  high! 


— JO«:. 
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The  time  histories  of  the  residuals  and  state  estimates  are  given 
in  Figure  13.  For  comparison,  the  residual  for  dor  obtained  by  setting  ■ 0 
is  given  in  Figure  14.  It  is  readily  apparent  that  the  scale  factor  error  could 
be  the  problem,  but  the  value  extracted  appears  unreasonably  high.  Of  course, 
the  scale  factor  could  be  accounting  for  other  unmodeled  error  sources  for  this 
record.  However,  this  does  support  early  suspicions  that  the  inconsistency  is 
indeed  in  the  measurements,  but  additional  analysis  would  be  required  to  ascer- 
tain the  exact  error  source. 


It  is  also  interesting  to  note  that  the  measurement  noise  statistics 
for  Ap  and  L9  in  Table  4 are  reasonable  values.  It  therefore  appears  that  the 
larger  values  obtained  from  the  R estimates  for  case  S-6  (refer  to  Table  8), 
could  possibly  be  due  to  the  assumption  of  white  noise  for  being  a poor  one. 

In  summary,  the  results  for  the  low  speed  records  correlate  favorably. 
It  is  clear,  however,  that  the  high  level  of  uncertainty  and  noise  present  in 
these  flight  records  precludes  the  use  of  visual  prediction  to  demonstrate  the 
accuracy  of  the  estimated  parameters,  that  is,  using  the  estimated  parameters 
from  one  record  to  predict  the  responses  of  the  airplane  to  another  record. 

It  is  believed  that  the  predicted  comparisons  would  be  highly  subjective  and 
a meaningful  criterion  by  which  to  interpret  the  results  hard  to  establish. 

For  example,  if  the  state  estimates  from  the  Kalman  filter  model  are  used  in 
such  a prediction,  the  high  level  of  process  noise  required  to  model  the 
control  input  uncertainty  would  clearly  "mask  out"  the  effects  in  the  predicted 
comparisons  introduced  by  the  differences  in  the  parameter  values  between  all 
three  flight  records. 

High  Speed  Records: 

Little  effort  was  devoted  to  the  high  speed  records.  As  expected, 
the  problems  cited  above  have  been  observed  to  be  considerably  more  serious 
for  these  records.  Structural  vibrations  appear  higher  for  these  records  and 
the  signal  to  noise  ratio  present  on  the  control  input  measurement  is  much  lower 
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than  the  low  speed  records,  although  the  characteristics  of  the  noise  appear 
similar.  Also,  lags  or  inconsistencies  in  the  instruments  take  on  added  impor- 
tance because  of  the  higher  bandwidth  of  the  airplane. 


From  the  knowledge  gained  with  the  low  speed  records,  it  was  decided 
to  drop  analysis  on  records  10  and  12  because  of  the  extremely  low  signal  level 
present  in  the  control  input  for  these  i >cords  and  identify  only  the  short  period 

j 

and  control  effectiveness  parameters  for  the  other  three  records.  Referring 
to  Table  2,  records  13  and  18  have  control  inputs  which  excite  principally 
the  short  period  motion  of  the  aircraft  whereas  record  19  has  persistent  pulse 
type  inputs  applied. 

The  identification  results  are  summarized  in  Table  9.  In  particular, 

I 

process  noise  was  used  to  model  the  control  input  noise  at  a level  selected  from 
the  PSD  plots.  Figure  1,  for  records  13  and  18,  and  the  R- estimator  used  to 
establish  the  level  of  measurement  noise.  The  &ax  and  (f  measurements  were  not 
used  for  reasons  cited  previously  and  because  the  phase  lag  in  ^ can  introduce 
large  errors.  For  record  19,  the  control  input  was  filtered  at  3.5  Hz,  process 
noise  was  assumed  equal  to  zero,  x^  identified  and  only  the  first  portion  of 
the  record  used.  Also,  R.  was  not  adjusted  to  be  consistent  with  the  R-esti- 
mator. 

The  time  histories  of  the  residuals  and  state  estimates  for  records 
13  and  18  are  shown  in  Figures  15  and  16,  respectively.  Response  comparisons 
with  the  flight  data  for  record  19  are  shown  in  Figure  17.  The  response  com- 
parisons appear  good.  However,  the  persistent  type  control  input  applied  for 
this  record  could  effectively  mask  out  or  hide  the  free  response  of  the  airplane. 

These  results  represent  the  very  first  efforts  to  analyze  the  high 
speed  flight  records  and  are  therefore  considered  preliminary,  at  best.  Time 
and  money  precluded  further  analysis.  However,  improvement  in  the  results  is  . 
very  unlikely  because  of  the  poor  quality  of  these  records. 


Section  VI 
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CONCLUSIONS  AND  RECOIN  DAT  IONS 

Flight  tost  data  from  a Navy  YT-21)  aircraft,  instrumented  and  operated 
by  NADC,  were  analyzed  to  estimate  the  longitudinal  stability  and  control 
parameters  of  the  YT-2B  at  two  flight  conditions:  a low  speed  landing  approach 

configuration  and  a higher  speed  cruise  configuration.  The  method  of  system 
identification  enployed  was  a nonlinear  iterated  Kalman  filter/fixed-point 
smoother  identification  algorithm.  The  method  also  had  the  capability  to 
adaptively  estimate  the  measurement  and  process  noise  level  from  the  flight 
data.  The  identification  results  obtained  were  good,  considering  the  extremely 
high  level  of  noise  in  the  flight  data. 

For  the  low  speed  flight  records,  the  parameter  estimates  correlated 
favorably  with  each  other  and  with  the  initial  wind  tunnel  estimates.  The 
covariance  matrix  of  the  parameter  estimates  also  appeared  to  be  a very  rea- 
sonable indication  of  the  estimation  accuracy.  Four  problems  were  responsible 
for  the  inaccuracies  in  the  estimated  parameters  and  the  differences  between 
the  parameters  extracted  from  different  flight  records.  They  are: 

1.  Inaccurate  knowledge  of  the  actual  control  input  applied. 

2.  An  inconsistency  between  the  angle  of  attack  and  normal 
stability  axis  linear  acceleration  measurements. 

3.  Incorrect  measurement,  caused  mainly  by  an  error  in 

the  trim  angle  of  attack. 

4.  Generally  very  high  level  of  noise  in  the  flight  records, 
due  to  the  flight  instrumentation  and  recording  system. 


Each  of  these  problems  was  discussed  in  the  main  text  of  the  report 
and  straightforward  solutions  posed.  The  inaccurate  measurement  of  the  control 
input  is  by  far  the  most  serious.  It  was  shown  that  both  filtering  of  the 
cantrol  input  and/or  modeling  the  error  in  the  input  with  process  noise,  where 
the  level  of  noise  is  established  from  the  sample  power  spectral  density  of  the 


24 


control  input  time  history,  were  very  effective  and  desirable  solutions.  The 
latter  solution,  of  course,  requires  that  the  identification  algorithm  be 
capable  of  accounting  for  both  response  error  (measurement  noise)  and  equation 
error  (process  noise).  It  is  clear,  however,  that  a properly  designed  and 
calibrated  instrument  at ion/ recording  system  would  eliminate  these  unnecessary 
efforts  at  the  analysis  level  and  improve  upon  the  accuracy  of  the  estimated 
parameters.  Because  of  the  generally  poor  signal  to  noise  ratio  of  the  control 
input  measurement,  meaningful  results  from  the  high  speed  records  were  not 
possible. 


The  following  -^commendations  can  be  drawn  based  upon  the  results  of 
this  investigation. 

1.  Data  analysis  and  flight  testing  should  be  done  concurrently. 
Problems  detected  from  the  data  analysis,  such  as,  modeling  or 
instrumentation  deficiencies  or  improper  control  input  excitation, 
for  example,  can  be  corrected  immediately.  To  do  otherwise,  can 
lead  to  a loss  of  valuable  flight  test  time. 

2.  The  instrumentation  problems  cited  above  should  be  investigated 
further;  initially  however,  at  the  flight  instrumentation  and 
recording  system  level.  Specific  flights  should  also  be 
devoted  to  ins .rumentation  calibration  checks.  If  necessary, 
the  combined  airplane/ instrument  at  ion  system  should  be  modeled. 

3.  A different  longitudinal  model  of  the  airplane  should  be  used 
for  further  data  analysis  at  these  flight  conditions.  The  model 
used  in  this  investigation  is  deficient  because  of  the  inherent 
inaccuracy  in  the  small  perturbation  body  to  stability  axis 
transformation,  especially  if  trim  angle  of  attack  is  not  known. 

4.  Further  analysis  should  be  done  to  obtain  estimates  for  the 
longitudinal  parameters  at  the  cruise  configuration  and  the 
lateral-directional  characteristics  at  low  angles  of  attack, 
but  only  after  the  instrumentation  problems  have  been  corrected. 


! 
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Figure  1 CONTROL  INPUT  POWER  SPECTRAL  DENSITIES  FOR  SELECTED  FLIGHT  RECORDS 


ACTUAL 


Figure  2 (cont'd)  ACTUAL  AND  FILTERED  CONTROL  INPUTS  FOR  RECORDS  2,  5 AND  6 


Figure  3a  RESPONSE  COMPARISONS  WITH  FLIGHT  DATA  - CASE  2-2 


Figure  3b  RESIDUALS  - CASE  2-2 


Figure  5a  (cont'd)  RESPONSE  COMPARISONS  WITH  FLIGHT  OATA  - CASE  5-2 


TIME-SEC 


Act 


TIME-SEC 


U Hz 


Figure  8 RESPONSE  COMPARISONS  WITH  FLIGHT  DATA  - CASE  6-2  SHORT  TIME  INTERVAL 


modified’  FC  - 2.0  Hz 


Figure  10  (cont'd)  RESPONSt  COMPARISONS  WITH  FLIGHT  DATA  - CASE  6-5,  Se  MODIFIED,  LONG  TIME  INTERVAL  (At 
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Figure  11  RESPONSE  COMPARISONS  WITH  FLIGHT  DATA  - CASE  6-5,  6 ACTUAL,  LONG  TIME  INTERVAL  (At  = 0.05) 


Figure  12  TIME  HISTORIES  OF  STATE  ESTIMATES  AND  RESIDUALS  - CASE  5-6  (Q  ^ 0 AND  <?e  UUFILTERED) 


CONSISTENCY  RELATIONSHIPS  - RECORD  5 


Figure  1 


TIME-5EC 
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(cont'd'  TIME  HISTORIES  OF  RESIDUALS  AND  STATE  ESTIMATES  FOR 
INSTRUMENTATION  CONSISTENCY  RELATIONSHIPS  - RECORD  5 
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Figure  14  RESIDUAL  FOR  Aa WITHOUT  SCALE  FACTOR  - RECORD  5 


TIME-SEC 


0 AND  <S„  UNFILTERED) 
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igure  17  (cont'd)  RESPONSE  COtdPARISONS  TO  FLIGHT  DATA 


TABLE  1 


MATHEMATICAL  MODELS 


Longitudinal  Equations  of  Motion: 
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TABU;  2 


SUMMARY  OF  FLIGHT  DATA 


Record 

Number 


tri  in 

dps) 


266.5 


trim 

(dot;) 


255.7  5.2 


251.6  5.9 


«V(.  Input  Type 


Oscillatory  type  bine 
wave  - short  period  only 

Oscillatory  type  - short 
period  only 

Pulse 


683.7 


669.9 


Pulse-short  period  only 


Pulse 


680.8 


682.8 


Double  doublet  type  - 
short  period  only 

Double  pulse  type  - short 
period  only 


Persistent  pulse  type 


* Pilot  indicated  moderate  turbulence  present 
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Para. 


Ms 

MS 


ze 

*2v 

*<x 


0. 

.0347 

-1.95 

-4.59 

-9.63 


0. 

■.0331 

0. 

.974 

.102 


TABLE  3 

INITIAL  PARAMETER  ESTIMATES 


/t  =236 

-.56147 

-.046 

-.1615 


Speed 

High  Speed 

=266.5 

EBST: 

Z&  =679 

-.56147 

-.56147 

-.56147 

-.0519 

-.049 

-.016 

-.205 

-.1723 

-1.2 

0. 

0. 

0. 

.031 

.032 

.113 

-2.20 

-2.09 

-4.96 

-5.85 

-5.30 

-18.6 

-12.28 

-11.11 

-51.8 

0. 

0. 

0. 

-.026 

-.0287 

-.004 

0. 

0. 

0. 

-1.242 

-1.124 

-2.46 

-.13 

-.118 

-.18 

* Parameters  assumed  known  sufficiently  accurate  to  fix 
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TABU.  4 


MEASUREMENT  NOISE  STATISTICS  FROM 
SAMIM.E  VARIANCE  COMPUTATION 


Standard  Deviations  (<J~ ) 


S 

.05  deg 

= 

1.4  fps 

= 

.1  dcg/sec 

°n 

= 

. 16  deg 

s 

3.46  ft/sec 

V 

= 

2.3  deg/sec 

= 

2 

.8  ft/sec 

= 

6.2  deg/sec' 

TABLE  5 

IDENTIFICATION  RESULTS  - RECORD  No 


Fix  - fixed  at  initial  estimate 


Fix  - fixed  at  initial  estimate  unless  otherwise  indicated 
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Appendix  A 

IDENTIFICATION  ALGORITHM 


This  appendix  summarizes  the  iterated  Kalman  fi ltcr/fixcd-point 
Smoother  Identification  Algorithm. 

Problem  Statement 

Consider  the  general  mathematical  model  representing  the  equations 
of  motion  of  the  aircraft  (or  other  dynamics  of  importance) 

X r f(x,  u)  -h  g(X,  f)  WiO, 

f = O (A-l) 

and  the  measurement  system  model  by 


y.  = h(xt  , f,ut)  + v't  , 


i=  Ot  1,  , N (A-2) 


where 


x - state  vector  for  the  aircraft  ( r/  , y ,0c.,  etc.) 
p - unknown  parameter  vector  ( M ^ * etc.) 

u - control  input  ( or  Xp  , etc.) 

- measurement  vector  of  sensor  outputs  at  discrete 
points  in  time 

f ( ) , h(  ) ( ) - functional  form  of  aircraft  and  measurement  system 

xo  - initial  condition 

and  *A/('t)  and  ixL  are  independent  random  noise  vectors,  which  arc  appropriately 
called  process  and  measurement  noise,  respectively,  and  are  chosen  as  white, 
Gaussian  and  zero  mean  for  analytical  simplicity  with  the  following  covariance 
matrices : 

E{i4f  (t  )*jvr(t+  T)]  = QS(T) 


£ (v'v'M 
III] 


= Kt 


U 
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Also,  the  uircraft  initial  condition,  x0  , and  initial  estimate  of  the  param- 
eter vector,  f0  , are  assumed  Gaussian  random  variables  with  covariance  matrix, 
Pa  ; which  is  a meusure  of  the  accuracy  of  the  "a  priori"  estimates  of  X0 
and  . Note  that  PQ  is  the  covariance  matrix  of  the  augmented  state  vector. 


The  problem  can  now  be  stated  formally:  given  the  measurement  data, 

, and  the  control  input,  ut  , for  i ■ 0,  1,  .....  N , obtain  the  best 
estimate  of  x0  and  , for  the  mathematical  model  chosen.  The  mathematical 
model  as  used  here  implies  both  the  structural  forms  for  f ( ),  q ( ),  and, 

h ( ) and  also  the  selected  values  for  the  <3  and  R.  covariance  matrices. 

The  best  estimate  will,  in  general,  depend  upon  the  criteria  used  to  determine 
the  best  estimate.  Bayesian  estimation  philosophy  is  used  here  and  the  best 
estimate  denotes  that  estimate  which  minimizes  the  mean  square  error  criterion 
(i.e.,  the  unbiased  and  minimum  variance  estimate). 

The  solution  to  the  problem  as  stated  above  is  that  of  obtaining  the 
fixed-point  smoothed  estimate  of  the  initial  augmented  state,  % = (xj  • P0T) 

which  is  the  conditional  expectation , £ where  Y(N)  = [yr,  y2r- •• 

With  appropriate  Gaussianity  assumptions,  this  problem  is  solved  by  using  an 
iterated  Kalman  filter/fixed-point  smoother  technique. 


Iterated  Kalman  Filter  Algorithm 


With  reference  to  Equations  (A-l)  and  (A-2) , the  iterated  Kalman 
filter  algorithm  is  summarized  below  between  the  arbitrary  data  points  at  the 
time  t-,}  and  t . Also  note  that  x is  interpreted  as  the  augmented  state 
vector,  i.e.,  x=  (xTPT)r. 


(j'O  /j)  ~ (j) 

s % , * y[>  • i . 

i/i  i / 1 - / t l 


A 


(augmented  state  estimates) 


x ■ , 

i/t  - / 


f (x,  u)dt  t- 


U) 


l.  L 


! -'A-' 


(predicted  state 
estimate) 


1 
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„<J)  . O')  .o)r  - 

p . , • r p . 4>  -f  Q 

* / 1 -I  i , t - / < - / 


(predicted  error 
covariance  matrix) 


,(j)  (j)  [ o)pM  fi(j)  1 (Kalinan  gain  matrix) 

t/i!  <■  | ' <' /<  1 > J 


A (}+l) 

XC-t/c 


\ ,,(j)  1 


X.  > P £ //  H P . //  * A?  • 2 

l-t/i-1  L-1  * L,  L-i  L [ <•  c/t-/  L I ^6 


(one-state  smoothed  estimate) 


pw  = rr  v'V'V,'1 

4 L " i J L/l  I 


(covariance  matrix  of 
state  estimation  errors) 


<j)  A Jf\ 
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At  + 


(transition  matrix) 


a t rorTAr 


C P Jt  , 


(linearized  output  matrix) 
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(predicted  residuals) 


A/  - sample  time  - t t:  / 


and  at  the  initial  or  starting  conditions 
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The  iteration  scheme  starts  with  j «1  and  terminates  when 

or  after  a prespecified  number  of  iterations,  usually  one  in  practice. 

■ <j)  A (J  O . „(J-0 


The  converged  values  of  x , 

A A ’ 


A 

mates  c 


•/* 

data  point. 


-Vi- 


and P. 


are  taken  as  the  esti- 


j,-  and  the  covariance  , respectively,  for  the  next 


Although  not  directly  a part  of  the  algorithm,  the  covariance  of 
the  predicted  residuals  is  given  by: 
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This  is  the  theoretically  calculated  covariance  of  the  residual  sequence  and 
it  can  be  used  to  detect  modelling  deficiencies  by  checking  to  sec  if  the 
actual  statistics  of  the  residuals  are  consistent  with  this  calculation. 


Fixed-Point  Smoother  Algorithm 
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Note  that  rQf-  is  not  needed  to  compute  xo,c  and  that  the  super- 
script (/)  denotes  the  last  value  of  the  particular  matrix  from  the  iterated 
Kalman  I liter  algorithm. 
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